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ABSTRACT 

We calculate the structure of accretion disks around Kerr black holes for ac- 
cretion rates M = 0.001 — 10M Q Such high-M disks are plausible candidates 
for the central engine of gamma-ray bursts. Our disk model is fully relativistic 
and treats accurately microphysics of the accreting matter: neutrino emissivity, 
opacity, electron degeneracy, and nuclear composition. The neutrino-cooled disk 
forms above a critical accretion rate M ign that depends on the black hole spin. 
The disk has the "ignition" radius r ign where neutrino flux rises dramatically, 
cooling becomes efficient, and the proton-to-nucleon ratio Y e drops. Other char- 
acteristic radii are r a where most of a-particles are disintegrated, r v where the 
disk becomes z/-opaque, and r tr where neutrinos get trapped and advected into 
the black hole. We find r a , r [gn , r u , r tT and show their dependence on M. We 
discuss the qualitative picture of accretion and present sample numerical models 
of the disk structure. All neutrino-cooled disks regulate themselves to a charac- 
teristic state such that: (1) electrons are mildly degenerate, (2) Y e ~ 0.1, and (3) 
neutrons dominate the pressure in the disk. 

Subject headings: accretion, accretion disks — dense matter — gamma rays: 
bursts 



1. INTRODUCTION 

A tight neutron-star binary loses orbital momentum via gravitational radiation and 
eventually merges, forming a black hole and a transient debris disk with a huge accretion 
rate, comparable to M & s _1 (see e.g. Ruffert et al. 1997 for numerical simulations). A similar 
accretion disk may form inside a collapsing massive star, so called "collapsar" (Woosley 1993; 
see McFadyen & Woosley 1999 for numerical simulations). In this case the central parts of 



1 Also at Astro-Space Center of Lebedev Physical Institute, Profsojuznaja 84/32, Moscow 117810, Russia 



- 2 - 



the stellar core quickly form a black hole of mass M ~ 2 — 3M that grows through accretion 
on the core-collapse timescale ~ 10 s. If the core material has specific angular momentum 
substantially above GM/c an accretion disk forms with M = 0.01 — 1M s _1 . 

Accretion disks are known to be efficient producers of relativistic jets in various sources 
associated with black holes, including X-ray binaries and quasars. It is reasonable to expect 
that the hyper-accreting disks in neutron-star mergers and collapsars produce relativistic 
jets as well. Because of the huge accretion rate, the power of such jets may be enormous. 
If a fraction ~ 10~ 3 of the accretion power Mc 2 is channelled to a jet, it will create an 
explosion with energy ~ 2 x lO 51 (M acc /M ) erg where M acc is the mass accreted through the 
disk. Hyper-accreting disk therefore provides a plausible mechanism for powerful relativistic 
explosions observed as gamma-ray bursts (GRBs, see Piran 2004 for a review). 

Accretion in the disk is driven by viscous stress t r( p that may be expressed as t r( p = ap 
where p is pressure inside the disk and a < 1 is a dimensionless parameter (Shakura & 
Sunyaev 1973). Viscous stress is created by magnetic fields that are amplified as a result 
of magneto-rotational instability. Numerical simulations of this instability show that a = 
0.01 — 0.1 (see Balbus & Hawley 1998 for a review). In this respect, the hyper-accreting disk 
is expected to be similar to normal accretion disks in X-ray binaries. It is, however, crucially 
different as regards microphysics of the accreting matter and its cooling. The optical depth 
to photon scattering is enormous and radiation cannot escape; it is advected by the matter 
flow into the black hole (or outward by the jet). The only possible cooling mechanism is 
neutrino emission. Significant neutrino losses occur at high M > 0.01M o s _1 , and then 
almost all the accretion energy is carried away by neutrinos. 

In any realistic scenario, the black hole has a significant spin as it forms from rotating 
matter and further spun up by accretion. The spin is likely to help the jet formation through, 
e.g. Blandford-Znajek process. It also increases the overall efficiency of accretion from 6% 
(zero spin a = 0) up to 42% (maximum spin a = 1). The black-hole spin has significant 
effects on the accretion disk because it changes dramatically the spacetime metric near the 
black hole, where most of accretion power is released. For example, in the extreme case of 
a — 1, the inner radius of the disk is reduced by a factor of 6 compared with the Schwarzschild 
case. This leads to a higher temperature and a much higher neutrino intensity. Therefore, 
disks around rapidly spinning black holes may create powerful jets via neutrino annihilation 
above the disk. 

The structure of neutrino-cooled disks was investigated in a number of works, and 
one of them studied accretion onto a spinning Kerr black hole (Popham, Woosley, & Fryer 
1999, hereafter PWF). PWF used accurate equations of relativistic hydrodynamics in Kerr 
spacetime, however, made simplifying assumptions about the state of accreting matter and 
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its neutrino emission, which later turned out invalid: (1) Electron degeneracy was neglected. 
The true degeneracy is significant: it strongly suppresses the e ± population, changes the 
equation of state and neutrino emissivity. (2) Neutron-to-proton ratio was not calculated 
and no distinction was made between neutrons and protons in the calculation of neutrino 
emissivity. The correct ratio is typically n n /n p ~ 10, which leads to a strong suppression of 
neutrino emission. (3) The produced neutrinos were assumed to escape freely. In the case 
of a Kerr black hole, this assumption breaks at M > 0.1M Q s" 1 — then the disk is opaque 
for neutrinos in the inner region where most of accretion power is released. 

Recent works discussed one or several of these issues. For example, Di Matteo, Perna, 
& Narayan (2002) focused on the effects of neutrino opacity at high M, however neglected 
degeneracy and assumed n n /n p = 1. Kohri & Mineshige (2002) studied strongly degenerate 
disks, and Kohri, Narayan, & Piran (2005) — the realistic midly degenerate disks. None 
of these works, however, attempted the construction of an accretion model in Kerr space- 
time — all assumed a non-rotating black hole and used a Newtonian or pseudo-Newtonian 
approximation. Pruet, Woosley, & Hoffman (2003) studied the Kerr disk models of PWF 
and pointed out that the disk must be neutron rich in the inner region, however did not 
attempt to model self-consistently the disk structure. Beloborodov (2003a) showed that f3- 
equilibrium is established in disks with M ^ lO 31 (a/O.l) 9 / 5 (M/M ) 6//5 g s _1 , where M is 
the black hole mass. The relation between temperature, density, and neutron richness in 
equilibrium was studied in that work, and the disk was shown to become neutron rich in the 
inner region if M > 10 31 (a/0.1)(M/M o ) 2 g s -1 . The structure of the disk was not, however, 
calculated. 

In the present paper we develop a self-consistent, fully relativistic model of accretion 
disks around Kerr black holes. Our study is limited to steady accretion with constant M, 
which is a good approximation at radii where the accretion timescale is shorter than the 
evolution timescale of the disk (it may be 1 — 10 s, depending on the concrete scenario of 
disk formation). The accretion rate is assumed to be constant with radius, M(r) = const, 
neglecting the fraction of M that may be lost to a jet. The model employs the custom- 
ary approximation of one-dimensional hydrodynamics (Shakura & Sunyaev 1973) where the 
effects of three-dimensional MHD turbulence are described by one viscosity parameter a. 

In § 2 we summarize the physics of neutrino-cooled disks, write down the set of relevant 
equations, and point out where progress is made compared with the previous works. We pay 
particular attention to the mechanism of cooling that couples to and regulates the state of 
accreting matter. The method of solution of the disk equations is described in § 3, and the 
results are presented in § 4. Our conclusions are summarized in § 5. 
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2. PHYSICS OF NEUTRINO-COOLED DISKS 

2.1. Outer Boundary Conditions 

Generally, the size of accretion disk is limited by the maximum angular momentum of 
the flow. We here consider only the region where accretion is quasi-steady, i.e. the accretion 
timescale is smaller than the evolution timescale of M. This sets an effective outer boundary 
r out which depends on the specific scenario of disk formation. In our models, we choose r out 
sufficiently far (well outside the neutrino ignition radius) to cover the whole neutrino-cooled 
region. 

Neutrino cooling is negligible outside the ignition radius r ign and the viscously dissipated 
energy is stored in the accretion flow, which makes the disk thick in the outer region. We 
choose the disk temperature at r out so that the energy content per baryon approximately 
equals virial energy ~ GMm p /r as is generally the case in any advective flow (e.g. Narayan 
& Yi 1994). The flow is assumed to be initially made of alpha particles at the outer boundary. 
The calculations will show that the flow completely forgets the boundary conditions as it 
approaches the ignition radius: a-particles are decomposed into free nucleons, and entropy 
per baryon and lepton number are regulated to certain values by neutrino emission. This 
convergence makes the details of advective accretion in the outer region unimportant. 

Note also that the thick disk in the outer advective region r ~ r out may produce a strong 
outflow, so that the accretion rate M decreases with radius. However, even this fact is not 
so important as long as we are given M at r > r ign where the disk becomes relatively thin 
and M remains constant. 

2.2. One-dimensional Relativistic Hydrodynamics 

The disk is described by vertically averaged quantities such as density p, temperature 
T, pressure P, energy density U, electron chemical potential /z e , neutrino chemical potential 
fj, u , etc. The disk is axially symmetric and steady, so all quantities depend of radius r only 
and we deal with a one-dimensional problem. The difference between radius in cylindrical 
and spherical coordinates is neglected as if the disk were geometrically thin. This difference 
is { r sph — r cyi)/ r ~ (l/2)(-^/ r ) 2 *C 1 where H(r) is half-thickness of the disk at radius r. The 
thin-disk approximation is quite accurate inside the ignition radius where H/r ^ 0.1 — 0.4, 
and less accurate at r ~ r out where H/r ~ 0.7 — 0.8. 

The vertically-averaged model is designed to describe most of the disk material, and 
it may not describe well the "skin" of the disk, especially if viscous heating is not uniform 
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in the vertical direction. The vertical structure of the disk is governed by the unknown 
vertical distribution of heating and the vertical neutrino transport (Sawyer 2003). We do 
not study the vertical structure in this paper and use the simple "one-zone" approximation 
of the vertically-averaged disk. 

The hydrodynamic equations of a relativistic disk express conservation of baryon num- 
ber, energy, and momentum (angular and radial) in Kerr spacetime (see Beloborodov 1999 
for a review). Hereafter we use Boyer-Lindquist coordinates x a = (t,r,8,<f)) with the cor- 
responding Kerr metric g a/ 3, which has two parameters: the black hole mass M and its 
dimensionless spin parameter < a < 1 (see e.g. Misner, Thorne, & Wheeler 1973). The 
disk is in the equatorial plane 9 = 7r/2 and the fluid motion is described by the four- velocity 
u a = dx a /dr = (w*, u r , 0, u^), where r is the proper time of the fluid. The equation of baryon 
conservation reads 

where the half thickness H is to be determined below by the equation of hydrostatic balance. 

The exact differential equations of azimuthal and radial motion may be replaced by much 
simpler and sufficiently accurate conditions following Shakura & Sunyaev (1973) and Page & 
Thorne (1974). This "thin-disk" approximation neglects terms ~ (H/r) 2 and assumes that 
the fluid is in Keplerian rotation with angular velocity, 



n = — = n K 



r 3 \ 1/2 GM 



+ a 



GM c 3 



(2) 



A small radial velocity is superimposed on this rotation, 



S \r 

where c s = [Pjp) 1 ! 2 is the isothermal sound speed and S(r) is a numerical factor. This 
factor is determined by the inner boundary condition (zero torque at the last stable orbit 
r ms ) and the Kerr metric; it is calculated in Appendix A. The accuracy of the thin-disk 
approximation is not perfect at large radii r ~ r out , where the disk is thick. However, the 
details of the outer region have no effect on the solution for the neutrino-cooled disk. 



The equation of energy conservation reads 

~ d{UH) (U + P) d(pH) 
dr p dr 



u 



(4) 



where F + and F are the rates of viscous heating and cooling per unit area of the disk, as 
measured by the local observer corotating with the disk. F~ depends on the state of the 
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disk matter that will be discussed in § 2.3 below. F + may be expressed in terms of Q and 
the kinematic viscosity coefficient u, 

F+ = vH(U + P) rg H {-g tt H (^) ' , (5) 

where 7 = (— g**)^ 1 / 2 ^* is the Lorentz factor of fluid measured in the frame of a local observer 
with zero angular momentum. The viscosity coefficient is related to the disk thickness H, 
sound speed c s , and a by v = (2/3)ac s H . 

The model described by equations (2), (3), and (4) neglects the effect of the stored 
energy and pressure on the radial and azimuthal dynamics of the disk, and retains the 
advected energy in the energy balance. This approximation is reasonable as may be seen 
from exact hydrodynamical models: the deviation from Keplerian rotation remains small 

10%) even when the advection effect dominates in the energy balance (Beloborodov 
1998). 1 Then the main effect of advection is the simple storage of energy that is not radiated 
away (F + — F~ 7^ in eq. 4). It strongly influences pressure and the scale-height of the disk, 
and hence changes u r according to equation (3). Equations (2) and (3) remain, however, 
good approximations. 



The vertical balance is given by 



Hy 2 JGMp 



r ) rP 



(6) 



where 

j = 2(r 2 -ar gv ^+0.75a 2 r 2 g ) 
2r 2 — 3r g r + ar g \/2r g r 
is the relativistic correction to the tidal force in the Kerr metric, and 

2GM 

r g = — . (8) 

In contrast to accretion disks in X-ray binaries and AGN, there is one more conservation 
law that must be taken into account in GRB disks. The lepton number (or, equivalently, the 
proton-to-baryon ratio Y e ) may change with radius because the neutrino and anti-neutrino 
fluxes from the disk may not be equal. The conservation of lepton number is expressed by 
the equation, 

(9) 



-^(iVp — N v ) — u r 



p dY e d , . 
+ -r{n u -n p ) 



m p dr dr 



1 A strong reduction of ft below fl K may occur in the limit of a large disk with no cooling. This limit 
does not apply to GRB disks. 
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where N v and N 9 are the number fluxes of neutrinos and anti-neutrinos per unit area (from 
one face of the disk), n u and n 9 are the number densities of neutrinos and anti- neutrinos 
inside the disk; they become significant when the disk is opaque (see § 2.4 below). Y e is 
related to the neutron-to-proton ratio by Y e = (n n /n p + 



2.3. Microphysics and Thermodynamic Quantities 



The disk is made of neutrons, protons, a-particles, electrons, positrons, photons, neu- 
trinos, and anti-neutrinos. The effect of magnetic field on the particle distribution functions 
may be neglected (see Beloborodov 2003a). The total pressure and energy density are given 
by 

P = P b + P 1 + P e - + P e+ +P V + P P , (10) 
U = U b + U 1 + U e -+U e+ + U u + U p . (11) 
Here, the baryon pressure and energy density are 



P b =-?-k B T(x f + - 



m. 



X 



U b = -P b , 



(12) 



where Xf is the mass fraction of free nucleons and 1 — Xf is the mass fraction of a-particles. 
Xf is found from the equation of nuclear statistical equilibrium (see e.g. Meyer 1994), 



4.9xl0 2 pro /2 T 1 i 9 o /4 exp 



16.4 



10 



4 



(1-*/ 



l-Y e - 



:i-x f r 1/2 , 



(13) 



where p w = p/10 10 g cm" 3 and T 10 = T/10 10 K. 

The radiation pressure and energy density are 

a T 4 
_ a r i _ 

1 7 3 ' ^' 
where a r = 7.56 x 10~ 15 erg cm " 3 K~ 4 is the radiation constant. 
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Electrons are neither non-degenerate nor strongly degenerate, and they are not ultra- 
relativistic at all radii. Therefore, no asymptotic expansions are valid, and the thermody- 
namics quantities for must be calculated using the exact Fermi-Dirac distribution. The 
pressure and energy density are given by the integrals, 



P,± = 



1 (m P c) 



3 TT 2 h 3 



-m P c 



r+oo 

Jo 



V 



dp, 



(15) 
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rr ( m eC) 3 2 

U e ± = — ^—m e c 



r+oo 

/ f(Vp 2 + hTVe)VP 2 + ±P 2 dp, (16) 
Jo 

where f(E,r)) is the Fermi-Dirac distribution function, 



f(E, V )= E ] (17) 
e e ' + 1 

9 = kT/m e c 2 , i] = fi/kT is the dimensionless degeneracy parameter, and ji is chemical 
potential. Since e ± are in equilibrium with radiation due to fast reactions e + + e~ ^7 + 7, 
and photons have zero chemical potential, one has the relation ji e - + /i e + = 0. We denote r\ e - 
by i] e and use r/ e + = —i] e in equations (10) and (11). The e ± population is then completely 
described by two parameters: 9 and i] e . 

The number densities of e~ and e + are 

(m e c) 3 r+co 



n P ± = 



7T 2 h 3 

The disk matter is neutral, which implies 

n e - - n e ^ 

This gives a relation between 9, r] e , p, and Y e 



r+00 

/ fWp 2 + l,TVe)p 2 dp. (18) 
^0 



= Y e ^. (19) 



The contribution of v and v to P and [/ becomes noticeable only in very opaque disks 
where neutrinos are completely thermalized and described by Fermi-Dirac distributions with 
chemical potentials p v and /ip = —p v (see § 2.4.2). f/^ and Up in the opaque disk are given by 
equations (27) and (30) below. The corresponding pressures are P u = U u /3 and Pp = Up/3. 



2.4. Cooling 

The cooling of the disk F~ may be written as a sum of three terms, 

F~ = F nuc + F U + Fp. (20) 

Here F nuc describes the consumption of heat by the disintegration of a particles as the flow 
approaches the black hole, 

dX 

F QUC = 6.8 x 10 28 p 10 — 1 u r H, (21) 
dr 

where all quantities are expressed in cgs units. This cooling is dominant in an extended 
region around 100r 9 where a-particles gradually disintegrate (see § 4 below). 
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The terms F v and F v represent the cooling due to emission of neutrinos and anti- 
neutrinos. There are four different channels of neutrino emission (see e.g. Kohri & Mineshige 
2002): (1) electron capture onto protons p + e~ — * n + v and positron capture onto neutrons 
n + e + — > p + 9, (2) pair annihilation e + + e~ — > i/ + P, (3) nucleon-nucleon bremsstrahlung 
n + n^n + n + v + v, and (4) plasmon decay 7 — > 1/ + v. The e 1 * 1 capture strongly domi- 
nates neutrino emission in GRB disks, and the other three channels may be safely neglected. 
The e ± annihilation into vv makes a small contribution even when electron degeneracy is 
neglected (PWF). When degeneracy is taken into account, the positron population is sup- 
pressed and the reaction rate becomes completely negligible. Bremsstrahlung and plasmon 
decay are important only at extremely high degeneracy and negligible in GRB disks. 

We calculate below F v and F p due to the e ± capture onto nucleons. This calculation is 
different in the transparent and opaque regions of the disk. 



2.4-1- Neutrino Emission from u -transparent Disk 



If the disk is transparent to the emitted neutrinos and anti-neutrinos, the emerging 
neutrino flux equals the vertically integrated emissivity that is found, e.g., in Shapiro & 
Teukolsky (1983). This gives, 

n r + °° 1 

E v = H Y e — K m e c 2 / f(E + Q,r) e )(E + Q) 2 [l - ————] 1/2 E 3 dE, (22) 



(E-QY 



1/2 



E 3 alE. (23) 



F p = H(l-Y e )-^Km e c 2 [ °° f(E - Q, -r] e )(E - Q) 2 
m P Jq+i 

Here Q = (m n — m p )/m e = 2.53, K = 6.5 x 10~ 4 s _1 , and f(E,rj) is the Fermi-Dirac 
distribution (eq. 17). 



Similarly, one finds the number flux of neutrinos, which will be used in the equation of 
lepton number conservation, 



m 



N v = HY e — K I f(E + Q, Ve )(E + Qf 



1 - 



1 



v Jo 



N D = H{l-Y e )-^K / f(E-Q,- Ve )(E-Qf 
m P Jq+i 



1 - 



(E + Q) 2 

1 



(E-Qy 



1/2 

E 2 dE, (24) 

n l/2 

E 2 dE. (25) 
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2.4-2. Neutrino Emission from v-opaque Disk 

Inside a zz-opaque disk, neutrinos relax to thermal equilibrium: a detailed balance is 
established between absorption and emission. Then neutrinos are described by the Fermi- 
Dirac distribution, and the energy flux of escaping neutrinos may be written as 

F - - T+k' (26 > 

where 

U v = ( -^rn e c 2 £°° f(E, Vu )E* dE (27) 

is the energy density of thermalized neutrinos inside the disk and t v is the total optical 
depth seen by u, including absorption and scattering (the cross sections of relevant processes 
are given in Appendix B). The chemical potential of neutrinos r\ v = p, u /kT that appears in 
equation (27) is related to /i e , /i p , and ji n because the detailed equilibrium v + n <-> e~ + p 
is established. The relation \i v + /i n = /i e + /i p then gives (see Beloborodov 2003a), 

^-^ = ln(i^+|. (28) 
When the disk is opaque also for anti-neutrinos, one finds 

F ° = ifl- (29) 

U, = J f(E, -^E* dE, (30) 

where we have used r\ v = —r\ v . 

Finally, the number fluxes of v and v that escape the disk are given in the opaque regime 

by 

n,..r. (m„r^ 3 r+ °° 



n v c (m e cY f + °° Nrf.p /oin 

(TT^)' ^ = wl f^EdE, (31) 

riyC {m e cf j" 

(1 + rp)' np ~2n^J 



N, = n p = §^ f°° f(E, - Vu )E 2 dE. (32) 



2.4-3. Transition Between Transparent and Opaque Regions 



Two changes happen at the transition to the opaque regime: (1) the probability of 
direct neutrino escape is reduced as (1 + t i/ )' 1 , and (2) neutrinos get thermalized. Note that 
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the rates of neutrino scattering and absorption are comparable, and a large optical depth 
t v ^> 1 implies that neutrinos are completely reabsorbed in the disk. The neutrino spectrum 
significantly changes at the transition: it is described by Fermi-Dirac distribution in the 
opaque region, with a non-zero chemical potential \i v . In particular, the mean energy of 
neutrinos, E u /kT, changes. 

We model the transition using the following approximate method. First we calculate 
the energy density of neutrinos that would be obtained in the transparent and opaque limits: 
jytransp _ p v j c f rom equation (22) and f/°P ac i ue from equation( 27), and define the parameter 



T rtransp 

x — v cyx\ 

^ytransp _|_ jy-opaquc ' V / 



If the disk is transparent we must find [/* ransp <C jj° p!Uinc (no reabsorption) and x — > 0. If 
it is opaque, [/* ransp 3> f/° paquc (strong reabsorption) and x — > 1. In our calculations, the 
neutrino flux is approximated by 

^transp c(l+7v) -l x < ±, 

[/ opaque c ( 1 + r ^-l j f x > 1 

This expression is continuous at the transition point x = 1/2. 

To model the change of N v at the transition we approximate the mean energy of the 
emitted neutrinos by 

ptransp popaque 

E v = (\-x)^- + x^ . (35) 

" v / ^y-transp jy-opaque \ / 

Then we define 

Nu = |K (36) 

which smoothly changes across the transition and has the correct limits at x = and x — 1. 
Transition to the opaque region for anti-neutrinos is calculated in the same way. 

Finally, we note that emission of muon and tau neutrinos is orders of magnitude weaker 
compared with the emission of electron neutrinos by e ± capture reactions. The main emission 
mechanism of and is nucleon-nucleon bremsstrahlung. The corresponding emissivity 
when nucleons are non-degenerate (p < 10 14 g s" 1 ) is given by (Thompson, Burrows, & 
Horvath 2000) 

^ - 10 30 CP? 4 (^) erg cm" 3 s" 1 , (37) 

where ( ~ 0.1 — 1 is a numerical factor. At high densities, when nucleons become degenerate, 
g^p M saturates at about 10 34 erg cm -3 s _1 . The muon-neutrino emissivity will be found to 
be well below F + / H and therefore may be neglected. 



2.5. Comparison with Previous Works 



Recent works on neutrino-cooled disks assumed a Schwarzschild black hole, and the 
main advantage of our model is that it is fully relativistic and describes accretion by Kerr 
black holes. Other advances compared with the most recent work by Kohri et al. (2005) are 
as follows. (1) Inclusion of the zero-torque boundary condition at the last stable orbit. This 
has a strong effect on the heating rate F + and the radial velocity u r . The effect of inner 
boundary is described by the factor S(r), which is much smaller than unity in the hottest 
region of the disk, S ~ 0.1. (2) Advection of heat and lepton number is treated accurately, 
using the differential equation of radial transport. This is essential since the disk is far from 
being self-similar and no analytical approximations are valid. The set of equations then 
becomes more complicated, however, it allows one to calculate the global model of the disk, 
from the outer advective region to the last stable orbit. (3) The transition to the opaque 
region is treated accurately. We find that there is only one thermalized neutrino species 
(electron neutrino), and its chemical potential fi u is obtained from the detailed equilibrium 
rather than assumed to be zero. (4) A corrected cross-section is used for anti-neutrino 
absorption by protons, which takes into account the proton recoil (see Appendix B). This 
minor refinement is interesting only in high-M disks that produce high-energy neutrinos. 



3. METHOD OF SOLUTION 

The disk is described by the set of coupled equations most of which are local, i.e. relate 
local parameters at a given radius. Two equations, however, are differential (eqs. 4 and 
9); they state conservation of energy and lepton number and contain advection terms that 
describe radial transport of energy and Y e . We therefore have to specify two boundary 
conditions. 

Our outer boundary is outside the neutrino-cooled disk, in the advective region, where 
neutrino emission may be neglected and matter is dominated by a-particles. So, one of our 
boundary conditions is Y e (r out ) = 0.5. 

As the other boundary condition one may specify any parameter that gives a reasonable 
approximation to the advection-dominated solution (Narayan, Piran, & Kumar 2001). For 
example, one could specify a certain value of (H/r) 2 ~ 1/2 that is characteristic for advective 
disks, and calculate all other parameters at r out using the local equations. There is a limited 
freedom in the choice of if(r out ), which reflects physical uncertainties in the behavior of 
the disk outside our boundary r out : how far the disk extends beyond r out and how much 
mass and energy it has lost to a wind. These uncertainties, however, have no impact on 
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the neutrino-cooled disk at r < r ign as we verify directly by varying the outer boundary 
condition: the solution we get at r ^ r ign <C r out is the same in all cases. Instead of H/r, one 
may use energy per baryon U/p as a free parameter at the outer boundary. In the advective 
disk, U/p is comparable to the virial specific energy GM/r. In the sample models shown 
below we specify U/p = GM/r at r = r out . 

The set of disk equations is a complicated mixture of differential and algebraic equa- 
tions, which involve integrals of the Fermi-Dirac distribution. We solve these equations on 
a logarithmic grid ri (i = 0..N), starting at tq = r out and moving inward. At each step, we 
need to find the parameters of the disk at r = V{ using the known parameters at rj_i > r, 
from the previous step. 

The state of matter at any radius is described by Y e and two independent thermodynamic 
quantities, which we choose to be 9 = kT/m e c 2 and i] e = p e /kT (because they enter as 
parameters in the integrals of Fermi-Dirac distribution). Density p, pressure P, and energy 
density U are expressed in terms of 9, i] e , and Y e as explained in § 2.3. For example, the 
expression for p is given by the charge neutrality equation (19). The computational problem 
now reduces to finding 9, i] e , and Y e at radius r\ given the known parameters at r^. 

Y e (ri) is easily found using the differential equation (9) since Y e (ri-i) is known from the 
previous step. The main difficulty is in the calculation of 0(rj) and r} e {rj). As two independent 
equations for 9 and rj e we choose hydrostatic balance (6) and energy equation (4). The scale- 
height H that appears in both equations may be expressed in terms of p using equations (1), 
(3), and (6), 

r?V = ?Af^r- (38) 



\r J a Anr 2 p \GMJ 
Then the hydrostatic equation becomes 



S GMM.A 1/3 
a Anr 3 J 



2/3 

(39) 



where r = r\. Both P and p are known functions of 9 and r) e , which are evaluated numerically. 
This gives one equation for 9 and r\ e . 



The second (energy) equation is differential and must be discretized, 



F + -F- = u r 



U t H t - U^H^ (U + P) ( Pt H t - p^HU, 
ri - p r-i- Ti_ x 



(40) 



All quantities taken at point r^-i are known and all quantities taken at are functions of 9 
and r) e . This gives the second equation. 
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We find 6 and r\ e that satisfy both equations (38) and (40) numerically, using a direct 
search in the (0,r) e ) plane. Once 9 and r\ e are found at r i; we calculate all parameters of the 
disk at this radius and move on to the next step r i+1 . 

4. RESULTS 

4.1. Sample Models 

Figures 1-8 show examples of the disk structure for accretion rate O.2M s _1 and three 
different values of viscosity parameter a = 0.1, 0.03, and 0.01. For comparison two cases 
are shown: a = (Schwarzschild) and a = 0.95 (Kerr). The black hole mass M = 3M is 
assumed in all models. 

First, we note significant differences between models with high a = 0.1 and low a = 0.01. 
Disks with low a accrete slower and have higher density. This has the following implications: 
(1) the region of neutrino emission in low-a disks extends to larger radii r ~ 200r 9 , (2) 
electron degeneracy rj e = \i e jkT is higher, and (3) Y e is lower. 

Electron degeneracy f] e is an important physical parameter that affects Y e , pressure, 
and neutrino cooling. In the outer advective region, r] e is decreasing as the heated fluid 
approaches the black hole until it reaches radii of a few hundred r g where nuclear cooling 
F nuc becomes significant. When only ~ 10% of a-particles are disintegrated (see Fig. 5), 
an energy Xf x 7.1 MeV^ 0.7 MeV is consumed per nucleon, which is comparable to the 
available heat stored in the flow at these radii. Then rj e begins to grow. 

The evolution of r\ e around 100r 9 is shaped by the competition of a few effects that are 
sensitive to a. In low-a disks, neutrino emission becomes significant at r ~ 200r 9 , which 
implies additional cooling (besides -F nuc ) and the drop of Y e from 0.5 toward a low equilibrium 
value. The coupled evolution of Y e , r] e , and neutrino emissivity reaches /^-equilibrium at 
r ~ 50r g . In high-a disks, neutrino emission is less efficient at r ~ 100r 9 and the drop of 
Y e occurs at smaller r, which leads to a different evolution of degeneracy r] e with radius. In 
all cases, however, /^-equilibrium is established with a low Y e as soon as neutrino emission 
becomes significant in the energy balance of the disk. 

The nuclear cooling becomes negligible where 90% of a-particles disappear (Fig. 5). In 
the a = 0.1 model, the neutrino cooling is still weak at this radius and the accreting fluid is 
quickly heated by F + . This recovery of disk heating leads to the knee in the profile of H/r 
at r w 40r g (Fig. 6). Then the growth of T above 1 MeV ignites strong neutrino emission 
that carries away the generated heat and produces the spike of F v {r-) between 10r 9 and 20r 9 , 
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1 10 100 1000 1 10 100 1000 

r/r g r/r g 

Fig. 1. — Temperature of the accretion disk with M = 0.2M Q s _1 around a black hole with 
mass M = 3M Q . Three disk models are shown with different viscosity parameters a = 0.1, 
0.03, and 0.01. Left panel: Schwarzschild black hole (a = 0). Right panel: Kerr black 
hole (a = 0.95). Radius is shown in units of r g = 2GM/c 2 , and temperature in units of 
m e c 2 = 0.511 MeV. 
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r/r g r/r g 
Fig. 2. — Density p(r) for the same disk models as in Fig. 1. 
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Fig. 4. — Y e (r) for the same disk models as in Fig. 1. 
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1 10 100 1000 1 10 100 1000 



r/r g r/r g 
Fig. 5. — Mass fraction of free nucleons Xf(r) for the same disk models as in Fig. 1. 
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Fig. 6. — Scale-height of the disk H(r)/r for the same disk models as in Fig. 1. 
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1000 



Fig. 7. — Ratio of neutrino flux F v + F D to the heating rate F + for the same disk models as 
in Fig. 1. 
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overshooting F + by a factor of 2 (Fig. 7). This overshooting effect does not happen in the 
low- a disks because they are cooled more efficiently around 100r r 

The Kerr disk models with a = 0.95 extend to a small radius r ms fa r g and have an 
extended neutrino-cooled region where F + xs F~ xs F u + F P) i.e. the approximate local 
balance between heating and cooling is established. This balance requires smaller scale- 
height H/r at smaller radii. As a result, H/r is reduced below 0.2. Significant differences 
between Schwarzschild and Kerr cases are observed in the inner region r < 10r g where most 
of accretion energy is released. The Kerr disk becomes opaque to neutrinos, and the value 
of Y e converges to Y e ~ 0.1 for all three viscosity parameters a = 0.1, 0.03, and 0.01. 

Figure 9 shows the contributions of P&, P 7 , P e = P e - +P e +, and P v to the total pressure 
P for two models of Kerr disk. One can see that the baryon pressure P& dominates in 
the neutrino-cooled region. This is a general property of all neutrino-cooled disks, which 
corresponds to their mild degeneracy i] e ~ 1 — 3. In the limit of strong electron degeneracy, 
P e would be dominant. In the limit of low degeneracy, P 7 ~ P e would be dominant. Only 
at mild degeneracy P 7 and P e are small compared with P b . The dominance of P b is thus 
a special feature of z/-cooled disks that is related to their self-regulation toward the mild 
degeneracy and low Y e ~ 0.1. 



4.2. Survey of Models 

A schematic picture of accretion disk is shown in Figure 10. The disk has 5 characteristic 
radii: 

1. - Radius r a where 50% of a-particles are decomposed into free nucleons. The 
destruction of a-particles consumes 7 MeV per nucleon, which makes the disk thinner. 

2. - "Ignition" radius r ign where neutrino emission switches on. At this point, the 
mean electron energy becomes comparable to (m n — m p )c 2 , enabling the capture reaction 
e~ + p — > n + v. Then neutrino cooling becomes significant, further reducing the disk 
thickness H/r. We choose the condition F v + F D = F + /2 as a formal definition of r; gn . 

3. — Radius r„ where the disk becomes opaque for neutrinos and they relax to a thermal 
distribution. Note that the disk is still almost transparent for anti-neutrinos at this radius. 

4. — Radius r p where the disk becomes opaque for anti-neutrinos, so that both v and 
v are now in thermal equilibrium with the matter. The disk is still cooled efficiently at this 
radius since v and V diffuse and escape the flow faster than it accretes into the black hole. 
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5. — Radius r tr where the timescale of neutrino diffusion out of the disk, t^m = (H/c)r u 
becomes longer than the accretion timescale, and neutrinos get trapped and advected into 
the black hole. The transition radius r tr is formally defined where F u + F p drops below F + /2. 

In addition, there is a radius beyond which the steady-disk model is inconsistent because 
of gravitational instability (e.g. Paczyhski 1978). We estimate this boundary from the 
condition 

^^ir 1 ' < 41 > 

The unstable region is where Q < 1. Besides the gravitational instability, the disk mass in 
this region becomes comparable to that of the black hole, and the gravitational potential is 
not described by the Kerr metric. 

We ran a series of models with various M, for two values of the spin parameter a = 
and a = 0.95. For each model, we found the five characteristic radii and the region of 
gravitational instability. The results are summarized in Figures 11 and 12. Contour plots of 
temperature T, density p, electron fraction Y e , and efficiency of local cooling (F u + F p )/F + 
are shown in Figures 13-16. 

As one can see from Figures 11 and 12, the radius of 50% disintegration of a-particles 
exists at all M of the sequence (M > 1O~ 3 M s _1 ) and weakly depends on M. In most 
models it is between 40 and 100r 9 . 

The ignition radius r ign exists if M > M ign which depends on a and a. Disks with 
M < M ign remain advective all the way to the black hole. For example, for the Kerr disk 
with a = 0.95 and a = 0.1, M ign O.O2M s _1 . By contrast, for the Schwarzschild disk 
with the same a = 0.1, M- lgn ps O.O7M s -1 . One can see in Figures 11 and 12 that the 
ignition radius first appears in the inner region when M = M ign . As M increases r ign moves 
to ~ 100r r 

The radii of transparency, r v and r p , scale with M approximately as M 3 / 2 , and r u Pi 2r p . 
Transparency of the disk for neutrinos also depends on the black-hole spin. The disk becomes 
opaque if M > M opaque ps 0.1M o s^ 1 in the Kerr case (a = 0.95) and if M > M opaque ps 
IMq s" 1 in the Schwarzschild case. 

Significant trapping of neutrinos occurs in the inner region of the Kerr disk if M > M trap 
which also depends on a and a. For example, for a — 0.1, M trap Pi 2M s -1 if a = 0.95 and 
M trap ^ 1OM s -1 if a = 0. The trapping radius r tr grows linearly with M. 

It is worth emphasizing that all three characteristic accretion rates, M ign , M opaquc , and 
Mt rap are lower for disks with smaller viscosity parameter a. The low-a disks are denser and 
have significantly larger r u and r tr . For example, the radius of opaqueness r v in the models 
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Fig. 9. — Contributions to total pressure P from baryons, electrons and positrons P e 



P e - + P e +, radiation P 7 , and neutrinos P v + Pp for the Kerr accretion disk 
M = O.2.M. b - 
with a = 0.01. 



= 0.95) with 

Left panel: model with viscosity parameter a = 0.1. Right panel: model 
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Fig. 10. — Schematic picture of the disk with characteristic radii indicated. 
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10 100 1 10 100 



r/ r g r/ r g 

Fig. 11. — Boundaries of different regions on the r-M plane for disks around a black hole 
of mass M = 3M and spin parameter a = (left) and 0.95 (right). Viscosity parameter 
a = 0.1 is assumed. Neutrino cooling is small in the shadowed region below the "z/-cooled" 
curve and above the "trapped" curve. The shadowed region marked "unstable" is excluded: 
the steady model is inconsistent in this region because of the gravitational instability. The 
disk extends down to the marginally stable orbit of radius r ms = 3r g for a = and r ms w r g 
for a = 0.95, where r g = 2GM/c 2 . 
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with a = 0.01 is almost one order of magnitude larger than in a = 0.1 models. 

The dependence of M ign , M op aquc, and M trap on a is well approximated by a power law 
(see Fig. 13), 

/ a \ 5 / 3 • / a \ • / a x 1 / 3 

-^ign = -^ign ^q - j^J ' -^opaque = -^opaque (^Q~yJ ' -^trap = -^trap \^g^J ' V^) 

where the normalization factors K depend on the black hole spin a. For a = 0, we find 
K lgIl = 0.071M Q s -1 , K opaquc = O.7M s _1 , and fsT trap = 9.3M Q s _1 . For a = 0.95, we find 
K ign = O.O21M s" 1 , ^opaque = O.O6M s" 1 , and K trap = 1.8M s" 1 . 

The maximum radiative efficiency of the disk, L u /Mc 2 , is determined by the binding 
energy at the last stable orbit r ms . It equals 0.057 for a = and 0.19 for a = 0.95. The real 
efficiency is somewhat smaller because part of the released energy remains stored in the disk 
and advected into the black hole. The efficiency is shown as a function of M in Figure 18. 



5. CONCLUSIONS 

This paper presents the model of neutrino-cooled accretion disks around rotating black 
holes that calculates self-consistently nuclear composition, neutrino emission, and fluid dy- 
namics in the Kerr metric. The model is one-dimensional in the sense that all parameters 
of the disk are integrated/averaged in the vertical direction and depend on r only, and the 
standard a-prescription is used for viscosity. The model also assumes a steady state, which 
is applicable only to radii where accretion timescale is shorter than the timescale of variation 
of M in the problem. The main advance of our work compared with PWF is the calculation 
of electron degeneracy and nuclear composition of the accreting matter; both dramatically 
affect the disk and its neutrino emission. 

The disk has a clear structure with five characteristic radii r a , r- lgn) r„, r p , and r tr . 
Radial advection of lepton number Y e and viscously dissipated heat is important outside the 
ignition radius r ign . Most of the viscous heat is lost to neutrino emission at r < r ign . 

The neutrino-cooled disk forms at accretion rates M > M ign which depends on the 
black-hole spin a and the viscosity parameter a (see Fig. 17 and eq. 42). General properties 
of the z/-cooled disk are as follows. 

1. — The disk is relatively thin, H/r ~ 0.1 — 0.3, especially in the inner region where 
most of accretion energy is released. The outer advective region r ^ 100r s is also significantly 
cooled by the gradual disintegration of a-particles, and its thickness is reduced. 

2. — The //-cooled disk is nearly in /^-equilibrium, in agreement with analytical estimates 
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(Beloborodov 2003a). In particular, the relation between p, T, and Y e calculated under the 
equilibrium assumption (Fig. 1 and 2 in Beloborodov 2003a) applies with a high accuracy. 

3. - - Degeneracy of electrons in the disk significantly suppresses the positron density 
n e +. However, the strong degeneracy limit is not applicable — the disk regulates itself to a 
mildly degenerate state. The reason of this regulation is the negative feedback of degeneracy 
on the cooling rate: higher degeneracy \i e jkT — > fewer electrons (lower Y e ) and positrons 
(n e +/n e - ~ e -Me/fer^ _^ wea k er neutrino emission — > lower cooling rate — > higher temperature 
— > lower degeneracy. 

4. — Pressure in //-cooled disks is dominated by baryons, P « P b = (p/m p )kT, most of 
which are neutrons. 

5. — All //-cooled disks are very neutron rich in the inner region, with Y e ~ 0.1 or lower. 

The high neutron richness has important implications for the global picture of GRB 
explosion (e.g. Derishev, Kocharovsky, & Kocharovsky 1999; Beloborodov 2003a,b; Rossi, 
Beloborodov, & Rees 2006). When the neutron-rich material is ejected in a relativistic jet, 
it develops a large Lorentz factor, and the neutrons gradually decay on scales up to 10 17 cm 
where the GRB blast wave is observed. 

The disks around rapidly rotating black holes are markedly different from disks around 
Schwarzschild black holes. They extend much closer to the center and reach higher tem- 
peratures and densities. For example, we found that the disk with a = 0.1 around a Kerr 
black hole with a = 0.95 becomes opaque for neutrinos at M opaquc ~ O.O7M s -1 , which is 
10 times lower than the corresponding M opaque for a Schwarzschild black hole. Besides, Kerr 
disks produce much higher neutrino fluxes in the inner region, with a higher mean energy 
per neutrino. The annihilation reaction v + v — > e + + deposits energy above the disk and 
can drive a powerful outflow that will be investigated elsewhere. 

This work was supported by NASA grant NAG5-13382. 



Conservation of energy and angular momentum is expressed by the equations (see Be- 
loborodov, Abramowicz, & Novikov 1997) 



Appendix A: Radial Velocity in a Relativistic Disk 




(43) 
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2 rut, 



where 



r ^ rr I 77 S 



(44) 



(45) 



is the shear and fx — (U + P)/pc 2 is the relativistic enthalpy per unit rest mass (// 1 in 
the neutrino-cooled disk); 7 = u*/ \/—g lt and f2 = -u^/m* are taken for the Keplerian circular 
motion (eq. 2). From equations (43) and (44) one can derive (see also Page & Thorne 1974) 
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MGM 



x 3 + a 



2tt c (x 3 -3x + 2a) 1 /%3/2 



3(xi — a) ; 



xi(xi - x 2 )(xi - x 3 ) 
3(x 3 - a) 2 



ln( 



In 



x 



Xi 



x 

, x ■ 



x 3 



3 

x Q ) - ^aln(— ) 

2 x 

3(x 2 - a) 2 
x 2 {x 2 - xi)(x 2 - x 3 ) "^x - x 2 



ln( 



x - x 2 



) 



x 3 (x 3 - x 1 )(x 3 - X 2 ) Xq-X 3 

where x = (rc 2 j 'GM) 1//2 , x± : x 2 ,x 3 are the three roots of equation x 3 — 3x + 2a = 0, and x 
corresponds to the marginally stable orbit r ms where fluid falls freely into the black hole and 
zero viscous torque is assumed. 



The radial velocity u r may now be expressed as 

M M va r ^ 

7T 



U 



2,rE x T <«> 
where cr^(r) and T(r) are known functions given by equations (45) and (46). Substituting 
the a prescription for the kinematic viscosity coefficient, v = (2/3)ac s H , one finds 



u' 



ac s 



fH\ 2Mra r 4> 



ac. 



(47) 



V r J 3tt T 

(37r /2)(T '/ 'Mra r .) varies from zero at the inner radius r ms to 



The numerical factor S(r) 
unity at r r g . 

"Newtonian" approximation that is often used in the literature on accretion disks (in- 
cluding GRB disks) is given by 

1/2 



Mr) = 1 " ( V 



It is derived for the accretion disk in Newtonian space by requiring conservation of angular 
momentum and imposing zero torque at a specified radius r ms , e.g. r ms = 3r g to mimic a 
Schwarzschild spacetime (Shakura & Sunyaev 1973). The correct function S differs signif- 
icantly from S N even for a Schwarzschild black hole: S ^ S N /2 in the inner region of the 
disk. 
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Appendix B: Cross Sections for Neutrino Interactions 

We summarize here the cross sections of neutrino reactions that we use in this paper 
(see Burrows & Thompson 2002 for a recent review of the reactions). The cross sections are 
expressed in terms of <t , 

= ^(mc 2 ) 2 c 1.71 x ID- 4 cm 2 . (48) 
7r(ftc) 4 v ; 

The neutrino energy is denoted by E and expressed in units of m e c 2 . 

1. — Neutrino absorption by nucleons: 

v + n — > e~ + p, z/ + p^e + + n (49) 

The cross section of z/ absorption by neutron is given by (e.g. Bemporad et al. 2002) 



'-<*w = j » + «y - (irfw • (50) 

where Q = (m n — m p )/m e = 2.53, #a — —1-26 is the axial coupling constant, (1 + 3g\)/A ~ 
1.44. This approximation neglects the recoil of neutron, however, it has only a small error 
< 1.5% when the neutrino energy is below 80 MeV (Strumia & Vissani 2003). For v ab- 
sorption by protons, however, there is a significant correction due to the recoil, and a better 
approximation should be used. We use the approximation of Strumia & Vissani (2003), 



a Pp (E p ) = 10" 43 n 2 (E p -Q j 2 J 1 



1 



k\E v - QY 
where k = 0.511. 



2. ■ - Neutrino-baryon elastic scattering. The cross sections of scattering on proton, 
neutron, and a particles are 



a p (E) = 



a E 2 



4 sin 4 9 W - 2 sin 2 9 W + 



0.30 a E 2 , 



<Tn{E) 



a E 2 fl + 3g 2 A 



0.36 a E 2 , 



a a (E) =4a sin 4 6 W E 2 = 0.21 CqE 2 , 



(51) 

(52) 
(53) 
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where 9w is the Weinberg angle and sin 2 ^ = 0.23. 

3. — Neutrino-electron (or neutrino-positron) scattering. The cross section of scattering 
for v or v is given by (Burrows & Thompson 2002), 



a e (E) = ^a 9E (l 



Ve 

4 



{Cy + C A f+ l -{Cy-C A f 



(54) 



where 9 = kT/(m e c 2 ) is temperature, r/ e = fx e /kT is the degeneracy parameter, C A = +1/2 
for v and C A = —1/2 for v e , and Cy = \ + 2 sin 2 9 W for both z/ and v. 

4. — Neutrino annihilation: u + u — > e~ +e + . Considering both the neutral and charged 
current reactions, the total cross section at high energies E V) E p 3> 1 is given by 



(T^p = X^cro — , (55) 

where P u and P p are the four-momenta of v and v in units of m e c, and = (1 + 4 sin 2 9yy + 
8sin 4 ^)/12 = 0.195 (Goodman et al. 1987; the full expression is found in Dicus 1972). 
After averaging over target distribution, one gets the average cross section for neutrino and 
anti-neutrino, 

a v (E v ) = ^K uP a E u E p (56) 

a- v {E v ) = ^K uP a E p E u (57) 
where E v and E 9 are the average energies of neutrinos and anti-neutrinos, respectively. 
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Fig. 12. — Same as Fig. 11 but for a = 0.01. 
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Fig. 13. — Contours of temperature T (in units of 10 10 K) on the r — M plane, (a) a = 0.01 
and a = 0. (b) a = 0.01 and a = 0.95. (c) a = 0.1 and a = 0. (d) a = 0.1 and a = 0.95. 
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Fig. 14. — Contours of density p (in units of g cm 3 ) on the r — M plane, (a) a = 0.01 and 
a = 0. (b) a = 0.01 and a = 0.95. (c) a = 0.1 and a = 0. (d) a = 0.1 and a = 0.95. 




Fig. 15. — Contours of Y e on the r — M plane, (a) a = 0.01 and a = 0. (b) a = 0.01 and 
a = 0.95. (c) a = 0.1 and a = 0. (d) a = 0.1 and a = 0.95. 
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Fig. 16. — Contours of local cooling efficiency (F u + F 9 )/F + on the r — M plane, (a) a = 0.01 
and a = 0. (b) a = 0.01 and a = 0.95. (c) a = 0.1 and a = 0. (d) a = 0.1 and a = 0.95. 
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Fig. 18. — Net efficiency of neutrino cooling of the disk L u /Mc 2 . The cases of a = 0.1 
and a = 0.01 are shown by the solid and dashed curves, respectively. The maximum effi- 
ciency that corresponds to complete cooling is shown by the horizontal dotted curve, (a) 
Schwarzschild disk (a = 0). (b) Kerr disk (a = 0.95). 



